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We study the R-parity violating minimal supergravity models accounting for the observed neutrino 
masses and mixing, which can be tested in future collider experiments. The bi-large mixing can 
be explained by allowing five dominant tri-linear couplings A' li2 ,3 and Ai,2. The desired ratio of 
the atmospheric and solar neutrino mass-squared differences can be obtained in a very limited 
parameter space where the tree-level contribution is tuned to be suppressed. In this allowed region, 
we quantify the correlation between the three neutrino mixing angles and the tri-linear R-parity 
violating couplings. Qualitatively, the relations |A^| < | A2 ] ~ | A3 1 , and |Ai| ~ | A2 1 are required by 
the large atmospheric neutrino mixing angle 82s and the small angle #13, and the large solar neutrino 
mixing angle #12, respectively. Such a prediction on the couplings can be tested in the next linear 
colliders by observing the branching ratios of the lightest supersymmetric particle (LSP). For the 
stau or the neutralino LSP, the ratio |Ai| 2 : | A2 1 2 : |Ai| 2 + | Aa | 2 can be measured by establishing 
Br{ev) : Br{[iv) : Br(ru) or Br(ue ± r^ 1 ) : Briy^r^) : Briy^r^ 1 ), respectively. The information 
on the couplings A^ can be drawn by measuring Br(litb) oc |A;| 2 if the neutralino LSP is heavier 
than the top quark. 

PACS numbers: 12.60.-i, 14.60.St 



I. INTRODUCTION 

The recent progress in neutrino experiments has established the picture for three active neutrino oscillations 0, 
which requires new physics beyond the Standard Model. The supsersymmetric standard model (SSM) would be the 
best motivated candidate for new physics. In the SSM, the lepton and baryon number conservation is not guaranteed 
by the gauge symmetry and thus one usually imposes an extra global symmetry, the R-parity, to forbid the fast proton 
decay. However, the proton longevity can be ensured by imposing only the baryon number symmetry. Then, the 
neutrino masses and mixing could be the consequence of the lepton number violation in the SSM. In this respect, the 
R-parity violating SSM has been extensively examined 0- 

The phenomenological study of the SSM strongly depends on the mechanism of supersymmetry (SUSY) breaking. 
The most popular scenario would have been the minimal supergravity model (mSUGRA) in which the universality 
of the soft supersymmetry breaking terms is introduced to solve the supersymmetric flavor problem. The mSUGRA 
model is highly predictive as it contains only 3 independent soft parameters: a universal gaugino mass Mi/ 2 , a universal 
scalar mass mo, a universal trilinear coupling Aq at the ultraviolet scale. The R-parity conserving mSUGRA model 
has been one of the most popular model as it contains a natural dark matter candidate of the Universe, the lightest 
supersymmetry particle (LSP), on top of its simplicity and predictability. Allowing R-parity violation in the mSUGRA 
model, such a feature is lost but the LSP can produce clean lepton number violating signals through its decay 
which provides a definite collider test of the model as the origin of the neutrino masses and mixing |J, yj 

This article aims at an extensive investigation of the mSUGRA model with R-parity violation. We will first analyze 
the parameter space of the model in which the phenomenological neutrino mass matrix consistent with the present 
neutrino experimental results can be accommodated. Concerning this, one may separate two issues: realizing (i) 
two large and one small mixing angles, and (ii) a mild hierarchy between the atmospheric and the solar neutrino 
mass scales. The property (i) can be explained simply by introducing an appropriate lepton flavor structure of the 
tri-linear R-parity violating couplings fA Il0|. Here, let us remark that the bi-large mixing can also be realized with 
allowing only the bilinear R-parity violating terms if one introduces non- universality in the soft terms |lll 1 1 121 . Il3f . It 
is however non-trivial to achieve (ii) as the mSUGRA model generically predicts very small values of Arn 2 o/ /Am 2 tm 
where the atmospheric mass splitting Am 2 tm comes basically from the tree-level contribution over-dominating the 
loop corrections which are supposedly responsible for the solar neutrino mass splitting Am^ ol . As a consequence, we 
find a very limited region of the parameter space consistent with neutrino experimental results. This is an undesirable 
feature of the mSUGRA scenario for the neutrino masses and mixing compared with other scenarios such as gauge- 
mediated SUSY breaking model 0. 

Given the allowed region of the (R-parity conserving and violating) parameter space, we will then examine various 
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ways to test the model in the future colliders through the decays of the LSP, which can be either a neutralino or a 
stau. As discussed in Refs. 0, HI El El j m models of neutrino masses and mixing with R-parity violation, one has 
specific predictions for various branching ratios of the LSP decay as the structure of lepton flavor violating couplings is 
dictated by the pattern of the neutrino mixing. For this, it is of course necessary that the LSP has a short lifetime to 
produce a bunch of decay signals inside the detector. Since the total LSP decay rate is proportional to the (heaviest) 
neutrino mass, the measurement of the LSP decay length could also be useful to test the model. In the parameter 
space accommodating the neutrino data, we find that the LSP decay length is shorter than a few centimeters. The 
neutrino mixing angles are basically determined by the ratios of the introduced trilinear R-parity violating couplings. 
We will discuss how strong correlations they have and how the ratios of the couplings can be determined by measuring 
the LSP branching ratios. As we will see, the quite reliable information of the neutrino mixing angles can be drawn 
from the collider experiments if tan (3 is not too large. For large tan /3, the close correlations between the mixing 
angles and the trilinear couplings are lost because of the large tau Yukawa coupling effect. 

This paper is organized as follows. In Sec. 2, we summarize how neutrino mass matrix can be generated up to 
1-loop level in the context of mSUGRA with generic R-parity violating interactions. In the analysis, we take the basis 
where the bilinear R-parity term is rotated away. In Sec. 3, we make a qualitative analysis to examine the sizes of 
various R-parity violating couplings in mSUGRA which are required to account for the current neutrino oscillation 
data. In Sec. 4, we calculate the LSP decay rate and branching ratios of various modes and find how the model can 
be tested in the collider experiments. We will conclude in Sec. 5. Various useful formulae for the 1-loop contributions 
to the neutrino masses and sneutrino vacuum expectation values will be given in the appendices. 

II. NEUTRINO MASS MATRIX AND FLAVOR STRUCTURE OF TRILINEAR COUPLINGS IN 

MSUGRA 

Let us begin by writing the superpotential in the basis where the bilinear term LiH 2 is rotated away : 

Wo = fiH r H 2 + hlUExEl + hjQiHxDl + h u k Q,H 2 U^ (1) 



W 1 = XiULzEl + X'iUQsDl (2) 

where Wq is R-parity conserving part and W\ is R-parity violating part. Here, we have taken only 5 trilinear couplings, 
A, and A^, assuming the usual hierarchy of Yukawa couplings. Among soft SUSY breaking terms, R-parity violating 
bilinear terms are given by 

V = m 2 L . Hl LiHl + BiLiH 2 + h.c, (3) 

where Bi is the dimension- two soft parameter. We will denote the Higgs bilinear term as BH\H 2 . In the mSUGRA 
model, the bilinear parameters, m\a 1 and Bi vanishes at the supersymmetry breaking mediation scale, and their 
non-zero values at the weak scale are generated through renormalization group evolution (RGE) which will be included 
in our numerical calculations. For the consistent calculation of the Higgs and slepton potential, we need to include 
the 1-loop contributions [ToL ITT| to the scalar potential as follows: 

*r l= 1 SteWln^-lf). (4) 



64tt 2 \ Q 2 2 

As is well-known, the electroweak symmetry breaking gives rise to a nontrivial vacuum expectation values of sneutrino 
(SVEV) Vi as follows: 

e __ (*«) _ rnl^+Bite + Yfcl 



where the 1-loop contributions E^ 1 ' 2 ' are given by E^. = dV\j H{dLi, Ej 2 ? = dV\ /L^dLj. Their explicit forms are 
presented in appendix B correcting various mistakes in the original calculation of Ref. 10]. The bilinear R-parity vio- 
lating parameters induce the mixing between the ordinary particles and superparticles, namely, neutrinos/neutralinos, 
charged leptons/charginos, neutral Higgs bosons/sneutrinos, as well as charged Higgs bosons/charged sleptons. The 
mixing between neutrinos and neutralinos particularly serves as the origin of the tree-level neutrino masses. We note 
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that the parameters £j should be very small to account for tiny neutrino masses. While the effect of such small 
parameters on the particle and sparticle mass spectra (apart from the neutrino sector) are negligible, they induce 
small but important R-parity violating vertices between the particles and sparticles, which in particular destabilizes 
the LSP together with the original trilinear couplings, Ai and A^. The derivation of the induced R-parity violating 
vertices has been presented in Ref. |7j. ^From the seesaw formulae associated with the heavy four neutralinos, we 
obtain the light tree-level neutrino mass matrix of the form ; 



M 2 

<r- = --^co S 2 /?, (6) 

where = MxM^jM^ + M§ cos 2(3/ fi with Ma = c^Mi + s 2 v M2- The R-parity violating vertices between particles 
and sparticles can give rise to 1-loop neutrino masses. Including all the 1-loop corrections, the loop mass matrix can 
be written as ^3 

M 2 

M t V = 'IT frSj+Si&coBP + Hij, (7) 

I'M 

where Uy denotes the 1-loop contribution of the neutrino self energy and 

/ -M 2 sm 2 6 w \ _ ( Mi cos 2 9 W \ 

yM-Mwt&nOw) M~M W J 

„ /sin/3\ „ (— cos/3\ 

+n ^{ -f J+n^(-^J. (8) 

Analytic expressions of the 1-loop contributions IT's are collected in Appendix A. Based on the neutrino mass matrix 
presented in the above, we will discuss whether the above mass matrix M u can account for both atmospheric and 
solar neutrino experimental data. 

As we will see, the tree level neutrino masses are much larger than the 1-loop contributions, which makes it difficult 
to account for the solar and atmospheric neutrino oscillations. As a result, the observed ratio Am 2 ol /Am 2 mt can be 
obtained in the mSUGRA only when some cancellation between the tree-level contributions to £, occurs to make the 
tree-level neutrino masses comparable to the 1-loop contributions or even smaller than the 1-loop contributions. 



III. NUMERICAL RESULTS: FITTING THE NEUTRINO DATA 



Let us now examine how the parameters in mSUGRA can be constrained by the recent results from neutrino 
experiments. For our numerical analysis, we scan the input parameters in the following ranges; 



lOOGeV < m < lOOOGeV, 
lOOGeV < M V 2 < lOOOGeV, 
OGeV < A < 700GeV, 
2 < tan/3 < 40. 



(9) 



The ranges for R— parity violating parameters we scan are 



4 x 1CT 6 < |Ai| < 6 x 1CT 4 , 
4 x 1CT 6 < |A 2 | < 6 x 1CT 4 , 



(10) 



3 x 10" 



4 x 10" 



4 x 10" 



< 



< 



< 



Ai 

a; 
a; 



< 10" 



< 10" 



< 10" 



We set the signs of Mi/ 2 and Aq arbitrary, but find that most of the allowed parameter space corresponds to the case 
that both signs are positive. 
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A. General behavior of neutrino oscillation parameters 

To understand the general feature of the R-parity violating mSUGRA model, we first plot the ratio A7n 2 oi /Am 2 tm 
for randomly generated points in the parameter space specified above. FIG. 1 shows that our model prefers the 
ratio in the range, 10~ 4 — 10~ 6 , which implies the dominance of the tree-level mass over the loop mass. Only a 
small fraction of the points are within the allowed range denoted by two horizontal lines. The big dots in this range 
are the solution points accounting for the bi-large mixing of the solar and atmospheric neutrino oscillations, which 
can be basically arranged by appropriate choice of the R-parity violating couplings. We note that more solutions 
can be found for tan/3 < 15. In FIG. 2, we plot sin 2 26*23 in terms of the ratios and |A 2 /A 3 | showing the 

correlations between the atmospheric neutrino mixing angle and the ratios of the R-parity violating parameters. The 
figure shows that the parameter^ has much cleaner correlation than A^ for the general points confirming the relation, 
sin 2 2023 = 4 1 ^2 1 2 1 1 2 / J2i l&| 2 M which is a consequence of the tree mass dominance. An interesting observation is 
that A^ also maintains a rather clean correlation, which implies that A^ gives the larger contribution to than A^. We 
also confirmed that the condition, |A^| < |A 2 | sa | A3 1 , is required to arrange the large atmospheric angle 6*23 and the 
small CHOOZ angle 6»i 3 . 

B. Behavior of the solution points 

In order to understand how the solution points can be obtained, we plot in FIG. 3, the quantity \m l ^ op — m 3 ree |/m 3 ree 
for the loop corrected and tree- level mass eigenvalue m 3 oop and m 3 ree , respectively. In the general parameter space 
(a), the ratio is shown to be smaller than 0.1. The situation becomes quite different for the solution points (b), for 
which the loop mass is comparable to, or even larger than, the tree mass. This implies that there must be some 
cancellation in & to reduce the tree mass. As a result, the clean correlation between and the angle #23 is lost, which 
makes it difficult to probe neutrino oscillation through the LSP decays, as will be discussed later. In fact, A^ will have 
the better correlation than ^ for the solution points. 

C. Importance of SVEV one-loop correction and neutral scalar loops 

In FIG. 4, we show how sizable the one-loop contributions to the SVEV are for the solution points. As can be seen 
from two figures, there occur almost order-one changes, which confirms cancellations in the tree-level values of £j. 
This implies that the loop correction to the SVEV should be taken properly into account for the determination of the 
neutrino oscillation parameters in mSUGRA models. 

It is also worthwhile to note that the the loop diagrams with virtual neutral scalars and neutralinos become sizable 
opposite to the case where the tree mass dominance works \j\ . In order to see this, we consider the parameters defined 
by 

* = (11) 

which are relevant for the neutral scalar loops 13] . In the general parameter space, one finds rji closely aligned with £j 
so that the neutral scalar loops do not give important contributions . In FIG . 5 , we plot | (^3 — 773) /?73 1 vs. | (£2 — ?72 ) /f?2 1 
for the solution points. The results indicate that rji can be much larger than £j and more importantly the alignment 
between £j and rji is lost in generic points (which are away from the diagonal line). All of these observations make it 
complicated to understand the structure of the neutrino mass matrix from mSUGRA models. Namely, all the loop 
corrections have to be properly included in the determination of the neutrino oscillation parameters. 

D. Fit to neutrino data 

Now we are ready to present how the trilinear R-parity violating couplings are constrained by the mixing angles 
and mass-squared differences observed in the atmospheric and solar neutrino oscillations. In FIG. 6 and FIG. 7, we 
show the correlations between the atmospheric neutrino mixing angle and the ratios £,2/ £,3 and A' 2 /A 3 , respectively. 
As indicated before, we find a better correlation with the latter. Even though A' 2 /A 3 has no analytic relation with #23 
for the solution points, we can obtain the following favorable ranges through the parameter scan: | A2/A3 1 from the 
neutrino data 

0.4 < IA2/A3I < 2.5 for tan/3 = 3-15 
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0.3 < \X' 2 /X' 3 \ < 3.3 for tan/3 = 30 - 40. (12) 

It is also amusing to find the correlation between the ratio A2/A3 and the solar neutrino mixing angle 612 as shown 
in FIG. 8. Similar to the case of the atmospheric neutrino oscillation, we can get the constraints: 

0.3 < |Ai/A 2 | < 1.6 for tan/3 = 3 - 15, 

0.2 < |Ai/A 2 | < 5.0 for tan/3 = 30 - 40. (13) 
In addition, fitting the measured mass-squared values, we find the allowed regions as follows: 

|Ai, a | , |A 2 , 3 | = (0.1-2) x HT 4 , (14) 

I a; 1 < 2.5 x io~ 5 . (15) 

Eqs. (12)-(15) are the key prediction of the mSUGRA model, some of which can be tested in the future colliders. 



IV. COLLIDER TEST OF THE MODEL 



In this section, we will investigate how the mSUGRA scenario for the neutrino masses and mixing can be tested 
in the collider experiments. Since the LSP is destabilized by the R-parity violating interactions, the structure of the 
R-parity violating couplings drawn in the above section may be probed by observing the lepton number violating 
signals of the LSP decay. Based on the parameter sets constrained by the neutrino data, we first determine the type 
of the LSP, which can be either a neutralino or a stau. Then, we calculate the cross section for the pair production of 
the LSP, and then its decay length and branching ratios. To consider the observability of the various decay modes, we 
will take the luminosity 1000/fb/yr in the next linear colliders [l4lll5| . For our presentation, we will show the trilinear 
R-parity violating parameters and the resulting neutrino oscillation parameters by selecting some typical examples 
of solution points and then give predictions of the production cross section and the decay length with important 
branching ratios of the LSP. As will be shown later, the LSP production cross sections are of the order 10-100 fb so 
that the branching ratios of the order 10~ 4 — 10 -3 will be measurable in the planned linear colliders. 



A. Stau LSP 



When the LSP is the stau, fi, it mainly decays into two leptons through the coupling A^. For small tan/3, t\ is 
almost the right-handed stau tr due to the small left-right mixing mass. Then the light stau almost decays into 
leptons via \LiL^E^ terms in the superpotential. Thus, one can expect that the branching ratios of those decay 
channels depend on the parameter A^. In this case, the following relation holds, 

Br{ev) : Brfjxv) : Br(ru) ~ | A x | 2 : |A 2 | 2 : |Ar| 2 + |A 2 | 2 , (16) 

as can be seen from Table I. Note that the decay length is much smaller than micro-meter (/im) so that the stau LSP 
production and decay occur instantaneously. The R-parity violating signals in the linear collider will be 

e + e~ — > fifx — > ifljvu 

, which are identical to the Standard Model background, 

e+e- -> W + W -> l+ljvv. 

But this is a flavor independent process and can be deduced from the total number of events to establish the flavor 
dependent quantities Br(liv). Therefore, the observation of the following relations: 

Br(rv) — Br(ev) + Br(fj,v), 
Br(ev) 



Br(fiu) 



0.1-2.6. (17) 



will be a strong indication of the R-parity violation predicting Eqs. (17) and (14). 

The situation is more complicated for large tan/3. A characteristic feature of this case is that there is a sizable 
fraction of the stau LSP decay into top and bottom quarks, if available kinematically, which is a consequence of a 
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large left-right stau mixing. However, the above relation l|16|l becomes obscured by the large tau Yukawa coupling 
effect. In TABLE II, we show this behavior for tan/3 = 38. One can see that the ratio of the branching ratios, 
Br(ev) : Br{^v) : Bt(tv) ~ 11.5 : 25.5 : 47.2, deviates from the \ coupling ratio, \\ Y \ 2 : |A 2 | 2 : |Ai| 2 + |A 2 | 2 ~ 15.6 : 
34.4 : 50.0. The deviation from the relation (17) comes from the large mixing between the stau and the charged Higgs. 
Another effect would be the charged Higgs contribution to the event: 

e + e~ -> H+H- -> t+t-vv . 



Even though a clean prediction for the r sector is lost, we are still able to establish the lepton number violating signals 
in the first two generations and measure the quantity: 



Br(ev) 
Br(fj,v) 



0.04- 25. 



(18) 



B. Neutralino LSP 



The lepton number violating signatures from the neutralino decay have been studied extensively as the LSP is a 
neutralino in the most parameter space 0,0.00- A characteristic feature of the neutralino LSP is that the vertex for 
the process Xi ~ > liW is proportional to which determines the tree-level neutrino mass (6). As a result, measuring 
the branching ratios Br(kjj) through either on-shell or off-shell W bosons will determine the ratio of |£i| 2 , that is, 

Br(ejj) : Br(fijj) : Br(rjj) = |£i| 2 : |f 2 | a : |& | 2 . (19) 

If the tree mass dominates over the loop mass, which is usually the case in the gauge-mediated supersymmetry 
breaking models 0, the neutrino mixing angles 623 and #13 can be cleanly measured in colliders Unfortunately, 
this is not the case in the mSUGRA model under consideration. As we discussed in the previous section (See Fig. 3), 
the tree mass has to be suppressed and thus the variables & are not correlated with the mixing angles 823 and #13 
in general and A^ maintain better correlations (See Fig. 6 and 7). ^From Tables III- VI, one can see that Br(lijj) or 
Br(liW) are proportional to £ 2 which, however, lost the correlation with the mixing angles. 

On the other hand, similarly to the stau LSP case, we can extract the information on Xi from the measurement of 
vlfr^ branching ratios for small tan/3 because the following relation, 

Briye*^) ■ Br{v^T^) : Br^r^) ~ |Ai| 2 : |A 2 | 2 : |Ai| 2 + |A 2 | 2 , (20) 

holds as shown in Tables III- VI. Likewise, if we measure the above branching ratios, the models can be tested by 
comparing them with the values allowed by the neutrino experimental data as in Eq. (18). Note that this is the case 
for small tan/3 < 15. For large tan/3, the relation (21) is invalidated again because of the large tau Yukawa coupling 
effect. 

Another interesting aspect is that the parameters AJ can be probed if the neutralino LSP is heavy enough to allow 
the decay modes, litb. The main contributions to these final states come from the couplings A^ and thus one obtains 
the following approximate relation: 

Br(efb) : Br(fitb) : Br(rtb) sa |A;| 2 : |A 2 | 2 : |A 3 | 2 (21) 

which should be consistent with the predictions (12), (14) and (15). In Table V, one finds the branching ratios for 
Utb large enough to be measured in the colliders with the integrated luminosity of 1000/fb. The branching ratios get 
too small if the LSP mass is more than twice top quark mass as shown in Table VI. This is because the vtt mode 
occurring through the light Higgs exchange becomes dominating. 



V. CONCLUSION 



We have investigated the neutrino masses and mixing originated from the R-parity violation and its possible collider 
test in the context of the minimal supergravity model (mSUGRA) where the soft terms are assumed to be universal 
at the ultraviolet scale. In this scheme, the bi-large mixing accounting for the atmospheric and solar neutrino data 
can be obtained by allowing five independent tri- linear couplings, A, and AJ. In the mSUGRA model, the tree- level 
contribution to the neutrino mass matrix is too much large compared with the one-loop contribution in the generic 
parameter space. In order to obtain the observed value of Am 2 oi /Am 2 tm , the parameters needs to be tuned to arrange 
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some cancellation in the tree-level mass which is made to be comparable to the loop mass. As a result, a nice relation 
between the R-parity (and lepton number) violating parameters and the neutrino oscillation parameters is lost. This 
is a dis-favorable feature of the mSUGRA model with R-parity violation as the origin of the observed neutrino mass 
matrix. However, in the allowed region of the parameter space, we still find the reasonable correlations between the 
tri-linear R-parity violating couplings and the neutrino mixing angles. That is, the bi-large mixing of the atmospheric 
and solar neutrino oscillation requires | | ~ | A3 1 (12) and |Ai| ~ |A 2 | (13), while the small mixing angle #13 requires 
|Ai|«|A^ 3 | (15). 

Such correlations predict some distinctive lepton flavor violating signatures of the LSP decay in the next linear 
colliders, which can provide a test for the R-parity violation scenario. If the stau is the LSP, there is no way to probe 
the couplings A^ and thus the two mixing angles #23 and #13. However, it is still possible to obtain the information 
on the solar neutrino mixing angle 6*i 2 since the stau decays into the final states through the couplings Aj. When 
tan/3 is not too large, the relation Br(ev) : Br{\iv) : Br(rv) — \\\\ 2 : |A 2 | 2 : |Ai| 2 + |A 2 | 2 can be checked to confirm or 
disprove our scenario. For large tan /3, the last relation in the r sector is obscured by the large tau Yukawa coupling 
effect. Thus, it will not be possible to establish the above relation without knowing other parameters like tan/3, 
Higgs masses, etc. For the case of the neutralino LSP, its decay into i/L can be observed to establish the similar 
relation as above: Br(ve ± T zf ) : i?r(^^ ± r =F ) : _Br(i^r ± T T ) = |Ai| 2 : |A 2 | 2 : |Ai| 2 + |A 2 | 2 . Again, such a clean relation 
is invalidated in the r sector for large tan (3. If the neutralino is heavier the top quark, it can decay into ktb through 
the couplings A-, predicting Br(etb) : Br(fitb) : Br(rtb) ~ lA'J 2 : |A 2 | 2 : |A 3 | 2 . Thus, its measurement can provide 
information on the atmospheric neutrino mixing angle. Let us note that such simple correlations are valid only for 
small tan/3 < 15. 
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APPENDIX A 



In the appendix A, we present the explicit forms of 1-loop contributions n's. Let us begin by fixing the notations for 
the diagonalization of the mass matrices. The mixing between the standard model particles and the supersymmctric 
ones in the R-parity violating model can be rotated away via the so-called seesaw diagonalization which makes the 
mass matrices block-diagonal, and then each blocks can separately be diagonalized. 



Neutrino- neutralino diagonalization 



The neutrinos/neutralinos mass matrix is written in 7 x 7 form as, 

M vN 
In M n 





(Vi- 




(3)- 







in the bases (u e , v^, v T , B, W3, H® , H®)- Rotating away the neutrino-neutralino mixing mass terms can be made by 
the following redefinition of neutrinos and neutralinos with 6 N 

nN~ \ 

(A.2) 

where (fj) and (x°) represent three neutrinos (v e , v^, u T ) and four neutralinos (B, W3, H®, if 2 ) in the flavor basis, 
respectively. The rotation parameters 8fj are given by 

N , M z ( s w M 2 c w M 1 M z M z 



{C > > ~ [c^M, + S * W M 2 ' " c^M, + s^M 2 ' ~ S0 ~' Cf) ~ ) ' ^ 

where Fn = MiM 2 /(c 2 v Mi + s 2 v M 2 ) + M|s 2/ g/^. Here sw — sm9w and Cw — cosOw with the weak mixing angle 
Ow ■ After rotating away the mixing terms, the mass matrix has block-diagonal form with 3x3 neutrino mass matrix 



8 



M v and 4x4 neutralino mass matrix in each blocks. In particular, the neutralino mass matrix is diagonalized by 
4x4 mixing matrix N, 

N^M N N = (M N ) dlag . (A.4) 



Charged lepton/chargino diagonalization 

Similar to the above case, charged leptons are also get mixed with charginos. Then the 5x5 matrix is 
written as 

in the bases (ej, W~ , Hi) , (e?, W + , ^ . Defining 9 L and 9 R as the two rotation matrices corresponding to the 
left-handed negatively and positively charged fermions, we have 



(A.6) 



where and denote the left-handed charged leptons and anti-leptons, (xj ) — ) an d (Xj) = (W /+ J ^2 )■ 

The rotation parameters d\- R are given by 

Oij = dcfcp - e t 5 j2 , Ofj = jr&cfcp and 
J F c (i 

(cf) = .^ m ^ te) .^f + A e) , (A . r) 

be A* [iMw M 

and i^c = M2 + My/S^pl fi. After block diagonalizing, the 2x2 chargino mass matrix M^±^± is diagonalized 

V^M^U = (M^ ± ) diag . (A.8) 



Sneutrino/neutral Higgs boson diagonalization 



Denoting the rotation matrix by 9 s , we get 



where 



S ttO 
il H l 



Hi 



y ;2 



(A.9) 



9 .s- 
9 s 

9f 3 



-& - ms^m\[m% z - m%.{m 2 A + M 2 z s 2 ) - m 2 A M 2 z s\c^\l 'Fs, 
+ViS /3 Ci 3 m 2 A [rr4 i - m\(m\ + M z cj,) + m A M z c 2 c 2 j3]/ 'Fs, 
-T]iS 2 p (? m 2 A M 2 z [ml. - m A c 2 p]/F Sl 
+r)iSpCpm 2 A M z [m 2 . + m A c 2 p]/Fs, 



(A.10) 



with Fs — {ni 2 . — mfj^m 2 . — m|^)(m?. — m A ) and m A ,mh and m# are the masses of pseudo- scalar, light and 
heavy neutral scalar Higgs bosons, respectively. Note that m 2 A = —Bn/cpsp in our convention. For our calculation, 
we assume that all the R-parity violating parameters are real and thus all 9's are real, too. We also note that the 
presence of the scalar fields as well as their complex conjugates in Eq. (|A.9(I is due to the electroweak symmetry 
breaking, which is expected to be suppressed with M z . Sncutrino and Higgs mass matrices can then separately be 
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diagonalized. 

Charged slepton/charged Higgs boson diagonalization 



Donning 9 C as the rotation matrix, we have 



~c* 



where 



0? 2 H 2 



of, = -6 



2 2 / 2 2 \ 



' (m 



AS 



9° 



111,-, 



s c & m A {m 2 Ri - m 2 H _) 



+Vi 



S l m A m Di 



( m H- - m e tl )i m H- ~ m U 

Sf}C m 2 A m 2 m 



+6 



- Zi 



l 2 H -)t/3 



( m H- - m2 eJ( m2 H- - m \ 2 Y 



(A.11) 



(A.12) 



( m H- - m lJ( m2 H - ~ m D ' wl ( m H- - m I 4 i)( m H- ~ m L) 

Here, m ff - stands for the charged-Higgs boson mass, and m 2 ^, m 2 ^ and m 2 Di correspond to the LL, RR and LR 
components of the i-th charged-slepton mass-squared matrix, respectively, and tn 2 n , 2 are their eigenvalues. We 

remark that the appearance of & in 8f x and 9f[ is due to the fact that we performed the rotation to remove the 
Goldstone mode from the redefined sncutrino fields. 



With the above notations, let us present the 1-loop contributions IPs: 



n« = 



^2{twN la -N 2 afm x oFf i 



gX 



+ 



lj3 



16tt 2 

gK x 



hiXj 

Si3 



16tt 

KK 

m 6 P 6£c 



0S m U la V ma m~-EO» - e^m T E0 T 3 
-eg n U la V ma m--mf + 0f 4 m r EO[ 
Og n U 2a V ma m.-Ql i - e^mrQli 
OfLU 2a V m am~- H0| - m r Qi 3 



16tt 2 

+ (i<— i) }, 



.9 



32tt 2 
5 2 iiy 



(A.13) 



jf^2^=[ ^ia^2am--K| i - 20&,E/i o V rno m x - HOf - <5 l3 2^ 4 m T HO- 



16tt 2 y^ 1 



U 2a V 2a m--d% i - 26? m U 2a V ma m -_COf 
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+ 0f m V ma U 2a m.-3 a - 5 a 9^ 2 m T y - 5 i3 6 j3 m r KT. } 

Kg ty, 



16tt 2 V2 
\ tw 
16tt 2 3\/2 



2m T COl - m r EO£ ] 
[2m 6 P 6o6o .+m 6 P 66 . ], (A.14) 



Q 2 

n jH?3 = - ^ (twN la - iV 2a )[ - N 3a m^oG^ + N^m^G^ } 
"' 1 =U la V 2a m--K^ 



16tt 2 V2 *• 

2, 

- S l3 03 2 m T J T - 5 i3 ^ 3 m r K5j ] 

=rr. 
2 

1 



16tt 2 V2 



™ 6 P 66 ,, (A.15) 



+ s i3 e^m T 3 T - s i3 e^m T mr 

- ^BfjJ 2a V ma m -_H0f 
+ ^[ fe032"vHOS + <5i 3 m r QI k ^ 3 ] 



' T^ 2 ""' 1 '^- (AJ(,) 



„2 



9 2 



16tt 2 * 
where the functions defined above are 

1 



j^U la V la m--K^ 

9 U 2a V la m.-Q^, (A. 17) 



y 2 * J 



BoK„,m^)-BoK„, m ^) 
+ZjiB {m\o,m 2 .) + ZijB {m\o,m\) 
+ ^e ih e jh B G {m\o,ml) + ^e iH e jH B {m\o,m 2 H ) 

~O iA e jA B (m 2 l o,rn 2 i ), 
G ti = 0gB o (mlo,m%) + ^s a 9 lh B o (mlo 7 m 2 ) 

-^c a 6 iH B (m 2 l o,m 2 H ) + ^s f:! e iA B (m 2 l o,m 2 A ) 7 
G 2i = °i2 s o(ra|o, ml. )- ic Q 6 ih B {m\o, m 2 h ) 

-^SaOiHBoim^Tn^) + ^s l3 e iA B (m 2 ,o,m 2 A ), 
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where 



Efj 



K5i 



QL 

p- - 

p- ~ * 



BO? 
BO[ 

HO? 
HO[ 
EOf 

Eor 



c 2 t B {m\- , ml n ) + S 2 5 (to 2 _ , mjj 

(egef lS 2 + egef 2 s 0C0 + egef^sp + egef 2 c}) B (m 2 , m 2 H - ) 



B (m 2 ,<)-*>« ,mlj 



+ (o?i0%4 + ege^spcp + ege%cp Sp + egef iC }) B (m 2 ^ , m \-\ 
{egc 2 - egasi) s (m|_ , mjj + (eg s 2 + egc lSl ) s ( m |_ , m\ a ) 



= {^A ~ Ogc lS% ) B (ml- , mjj + (6gs 2 + 9gc lSl ) B (m 2 ^ , m? J 

- (0«s/jc,j + #gc|) So (to 2 - , ), 
= (-flgciSi + egs 2 ) B (m 2 ^ , to 2 J + (Ogsid + 6g 3 c 2 ) B (m 2 ^ , to 2 J 

- (tias} + egspcp) B (m| _ , to^_ ), 
= Qii( m L- m r) ; 

= (-flgciSi + egs 2 ) B (m|_ , to 2 J + (flgs iCi + ^c 2 ) B (to 2 _ , to 2 J 

- + ^c 2 ) B (to 2 _,to^_), 
= -c b s b B (ml,m 2 J + s b c b B (ml,m 2 2 ) 

slBoim^mlJ + c 2 b B (m 2 b ,ml 2 ) 

„2n / 2 2 \ i 2 o / 2 2 \ 



c b B (m b ,m- bi ) + s b B (m b ,m~ b2j 



s 2 B (m 2 , m\ n ) + c 2 B (m 2 , mjj 



+ + og 3 ef 4 s 0C/3 + eg 4 e%cp S0 + egef 4 c 2 ) s (m 2 _ , m 2 H 

s 2 B o (m 2 ,m 2 H -) 



+6gef 1 S ij c 2 B (m 2 , m? a ) + s 2 B (m 2 , m? i2 ; 



So(m 2 ,TO? sl )-i3o(m 2 ,to| 2 ) 



c 2 £ (m 2 , mf^ ) + s 2 B (m 2 , m? i2 ) 



B0?(to 2 

i \ x 



m T ), 

^(<-.<)-^K.-.<) 



i 2 - 

Xa 



HO?(m 

c i B o(m^-,™L) + s i B o(m~-,m| i2 ) 
EOf(m|_ - m T ), 



ft/, - 



0*1 



-*a(fifi + ^) + Ca (^ + ^; 



2 TO-.C ct _^-Af z C 2/ 3C Q+/ 3 



c . 2 m I, S 



("*!, - m 2 ) (ml. - m 2 H Y 
— M|c 2/ 3S Q+ ^ 



(m|. -to 2 )(to? s -mfy)' 



TO, 



(A.18) 
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Zij = rjiiijm^M^CpSp 



ml. 



n + n 



(A.19) 



Here the angle a is the usual diagonalization angle of two CP even Higgs bosons, and F s = (rn 2 i —m 2 i )(rnl.—'m 2 i )(rn 2 , i — 
m 2 H ). Recall that the angle a is defined by c 2a — C2p{m 2 A — M§)/(m 2 — m 2 H ) and s 2q = S2p{m 2 A + M§)/(m 2 — m 2 H ). 
Ci = cos 9i, Si — sin 6i are the components of the diagonalization matrices of the slepton mass square matrices. 
i = e, /i,f, each. c&,S(, are the same ones for sbottom mass square matrices. B is the loop function defined as 
Bo(z,y) = -^lnf-ln^ + l. 

APPENDIX B 

We collect here 1-loop contributions to the sneutrino minimization condition by calculating all the field-dependent 
particle masses. As given in the paragraph, the sneutrino vacuum expectation values are given, 

, _ (Oi) ml Hi +B t t + ^ 

where the 1-loop contributions T, L ! are given by e£? = SVi/HldLi, s£ } = dV 1 /L*dL i . Here, we will give the 
explicit forms of these Sl/s. 

1. Top squark contribution 

The contributions from the stop mass eigenstates as follows: 



- ^ o 2 4M h 2 T m o K f m ^ - 1 ) > (B.2) 

LiK3> 64tt 2 1 2m? - m? - m? v 2 s 2/3 /2 J Q 2 1 

y ( 2 ) rf ^_ 6 [ 1 M 2 M 2 Lt - M 2 m 8M 2 , - 5M 2 I / ™g, 

S ^ fe) - ^ i 2~ + 2m? -m? - m? 6^ f ( "Q 2 " ~ 1 ) ' (R3) 



,2 



where M 2 t = M 2 n , M|, t = -M~ 12 and M^ t = Al 2 22 ignoring small R-parity violating contributions. 

2. Bottom squark contribution 



The sbottom contributions are 



647r 2 I 2 vi 2m 2 — m? — m~ v\ 

bi 



2 



where M£ fe , A/j^, and Mp b are defined as in the stop case. The contributions from the first two squark generations 
can be obtained by obvious substitutions. 

3. Charged slepton or charged Higgs contributions 

Let us denote the charged slepton and Higgs contributions to E's as 
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where 4> stands for the mass eigenstates en,ei2 and H . Here S^\if>) = (dm^/diij) /2v\ and S^\(j)) = 
(dm^ ) /duj)/2uj. M 2 . = M 2 tll ,Mp. = M 2 il2 and M^. = M 2 22 . Where i runs e,fl,f. ignoring small R-parity 
violating contributions. 



For i = 1,2, 



For i = 3, 



For i = 1,2, 



For i = 3, 



4. Charged Higgs contributions 



v 2 (Ml t -m 2 H _) 



o(l)^-^ _ (c^3"W, gl )(K- -M| 3 )(m^-2M^)s^ + /xm T M|) 3 ) 



5 * {H ) = W{m* H _-Ml) ■ (R9) 



si 2 \H-) 

= { ((m 2 ff _ - m^fm^ - m^M 2 , - M 2 - m 2 )c 2/3 

- 2(m 2 H _ - M 2 R3 ) S 2 2f) {M^ + M^m 2 + ra\) - 2(m 2 H _ - M 2 3 )c 2 ^ 2 m 2 T 

+ 4 S/3 c^ M m T (-m 2 + 2M^)MlJ / 2v 2 (m 2 H - - m 2 fl )(m 2 H - - m? 2 ). (B.10) 

5. Scalar Lepton contributions 



1 (6L) " v 2 M 2 Li {m 2 H _-M 2 Li ) ' (BJ1) 

S${e L )=0. (B.12) 

5{ 1) (e fl ) = 0, i = 1,3 (B.13) 

v M^(^(m^_C2^ + M 2 2 )i? 2 + (m 2 ^ - M 2 Jm 2 2gi ) 
* 2 ^ = v m 2 2 (m 2 H _ - Ml) ' 

sS(Al) = 0. (B.15) 

S$ L) (fi R ) = 0, * = 1,3 (B.16) 



c (2),- , _ «- - M 2 Li )(M 2 M 2 Li + 2M^) - 2M^_^ 
1 CL '~ 2v*M 2 Li {m 2 H _ -M 2 Li ) ' [ ^ U) 
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For i = 1,2 



For i = 3, 



where, 



For i = 1,2, 



sgfe) = Z 2M|+M| (B18) 
S f'(«) = il#, i = 1,3 (B.19) 



S<>>(h, 2 ) = -— A.T-^rA, (B.23) 

f 777- — 777- 

T2 Tl 



S3 (^1,2) = r-5 o 2~7 — 2 TTi — 1 2 \ ' (B.24) 



(A) = (m 2 H - - m\ ) (m 2 ^ - m\ ) (M 2 3 ,im T + M\ 3 spM^) 

+ m%- ( M k( M k - m H- ) + M£ 3 )(2M& - cpm 2 T )spc 
+ M D s m H- ( m f 2 + - m 2 H -)fim T c 2 p, 

(B) = (M 2 R3 m 4 H - - i<M! 3 + M^M*, 

- ^4 3 m^_ + M 2 3 M 2 R3 m 2 H - - 2M% 3 m 2 i -)c fi M 2 , 

- m \ m \ 04- - M R 3 ) c P m r + M D 3 m H- ( m H- - m \ - ml^spcpfimr 
+ ™V (M 2 R3 (M 2 R3 -m 2 H -) + M£,)(2M& - m 2 T )c% 

(C) = - (m^_ - m 2 fl ){m 2 H - - mf 2 ) 

x(M 2 D3 (m 2 fi + m\)pLm T + (2M% 3 + M 2 3 (M 2 3 - M| 3 ))^M^) 

- Ml 3 m 2 H - {{m 2 fl + m? 2 )(m? 1 + m? 2 - ) - 2m\m\ 2 )(?^m r 
+ ™ 2 h- (M e R3 (M 2 3 + m 2 H - )M Rs + (m 2 H - M 2 3 + 3M 4 D 3 )M 2 R3 

+ (M 2 3 2 m 2 H ^)M i D3 )(m 2 T - 2M 2 v )spc% 
(£>)= - m? i m? 2 ((M2 3 -Mi 3 )(M2 3 -m 2 ff _)+2M| )3 )c^ 

+ (-2M* 3 + (-M^ + AM 2 R m 2 H _ - 2m^ + ZM 2 L M 2 R3 )M A D3 

+M 2 R3 (M 2 3 + m 2 H -)(M 2 R3 - m 2 H .)(M 2 R3 - M 2 3 ))c fj M^ 

+ Mn 3 m 2 H -({m 2 fl + 777| 2 )(777? 1 + 777? 2 ~ m\-) - 2777^ 777^ )spCf3[im T 

+ ( Ml 3 - (Mf 3 + m 2 ^ )M* 3 + (tt7 2 h - Af 2 3 + 3M 4 D3 )M 2 R3 

+ (M 2 3 - 2ra 2 H - )M% 3 ){m 2 T - M 2 ^ . (B.25) 



S i N-^(1T MfK -ml) J' (R26) 
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For i 



where. 



4 2) (ri,2) 



(F) 



T 2 



IT 



Av 2 m% m% (m 2 _ — m% )(m 2 _ — m% ) ' 



(B.27) 



(E) = 



(F) 



(in 



2 \l 2 2 \ 

m f 1 ){ m H- - m f 2 ) 



M 



x {m\m\ (2M 2 , - M|) + 2M\^m\ + AM^M^nirSp - 2M 2 R3 M^) 



AM 



d. 



.m 2 H - (m fl + m? 2 - m 2 H ^){jiml - 2 y um T M^)s /3 c| 



8m 2 H - {M% 3 + M Rs - M 2 R3 m 2 H -)(c 2 s 2 M^ - c%m\ + c%M 2 w ml) 
2m 2 H - (M 4 D3 + Mi 3 - M\ 3 m\- )c 2 ^ 2 m 2 T 

4((2m^_ - M 2 3 )M A D3 + M 2 3 (M 2 3 + m 2 H .)(M 2 R3 - m 2 H _ ))c 2 f} M 2 f m 2 T 
2m 2 fi m 2 f2 (M 2 3 -m 2 H -)c 2 m i T , 

( m H- - m f 1 )( m H- - m \) 

m\m\{M 2 3 - M 2 3 )(iM 2 - 2M 2 N ) 

-2{m\m\- M* D3 - M^m 2 . 
+4M 2 D3 (M 2 3 + M^s^M^mr 

-2(m 2 fi mj 2 -M 4 D3 ^M R3 )M^' 

■ TO f 2 X m fi + m \ - m H- ) - 2m \ m X) 

x (M 2 3 M 2 R3 - m 4 H - )c 2 m 4 T - 2M 2 3 m 2 H .Sfsc 2 0f i(m 3 r - 2M 2 v m T ) 
2m 2 H - {s 2 fj Ml + c^AM^ml - m*)) 



2((m? 



fi 3 



(ML 



A- )M R3 + (m 2 H - M 2 3 + ZM% 3 )M 2 R3 + M 4 Di (M 2 3 2m\- ) 



o.™2 2 2 2 

2m H -Cg/j, m T 



Mr 



« + m%- )Mi 3 + (m|_ M 2 R3 + ZM^M 2 ^ + M 4 Di (M 2 Ra - 2m 2 H - ) 



+ 



8M« s + 4(M* 3 - (3M| 3 + 4m 2 ,- )M 2 3 + 2m 4 ff _ )M* 3 

-AM 2 R 3 (^3 + "4" ) ( M «3 - ) ( M «3 - ™H- ) 



(B.28) 



6. Neutrino-Neutralino contributions 



The sneutrino VEVs contributions arc 



(B.29) 



where -0 runs for four neutralinos x°. For 4 ncutralinos, if/ = = 1,4, 



2^M| 



- nM^m% + {-s 2f) M 2 z M^ + ii(Ml<? w + M%s 2 w ))m% 
+ M^(s 2 /3M|M^ + /i(M|+/i))m* 

- (M 2 M? + ^ 2 (M 2 c 2 w + M 2 s 2 w )) 1 / £>(*), 



(B.30) 



where, 
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8. CP-odd sneutrino or Higgs boson contributions 
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{<t>) = 



,(1,2) 



64tt 2 



{<t>)ml ln-^-1 



Q 2 



where <f> runs for Tm(v), A . Here, are given by 

Sj 1) (i?/)=0 ) 



^(A°)=0, 

5f)(A0) = -I^cos 2/3 . 



(B.40) 



(B.41) 
(B.42) 
(B.43) 
(B.44) 



9. Charged lepton or chargino contribution 

Neglecting the contributions from the two light charged leptons by taking m e /J = 0, we get 



^(*) = -^ (1 ' a) (*)^(^^-l)- 
where \& runs for r and the charginos x]~ 2 an d S^ 1 ' 2 ^^) are given by 



(B.45) 



S , f 1) (*) = -A i m T ^ (m|-m^ + )K-mL) /£>(*) 



Vl 



X2 



(B.46) 



where, 



4 v 



+ mJ) + mX(/i 2 + 2M^4) 

+ %m?{ mV2 C 2 + m|[M^(3 - t%) - Mf /<$] 

- 2M% r {M 2 ntp - 3M^s} - v 2 )} } /£>(*), 



£)(*) = 5m| - 4m^,(m 2 _ + m 2 _ + m 2 ) 

Xi X2 



Xi 



x 2 



Xi X2 



Xi X2 



(B.47) 



(B.48) 
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FIG. 2: The atmospheric neutrino mixing angle vs. and |A 2 /A 3 | for general points. 
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FIG. 3: The ratio \rri-j, — m%\/m% for the loop-corrected (7713) and tree-level (m§) neutrino masses with varying A2. 
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FIG. 4: Plots for (£3 — £,z)/£,3 versus (£2 — £2)/ £2 showing the effect of SVEV 1-loop correction for solution points 
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FIG. 5: Plots for |(^3 — 773 ) / ^3 1 versus |(^2 — 772) /^2 j showing the deviation of rji from £» direction for solution points. 
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(a) sin 2 2e 23 vs. |£ 2 /£ 3 | 



(b) sin 2 20 2 3 vs. A 2 /A 3 



FIG. 6: The correlation between the atmospheric neutrino mixing angle and ^2/^3! or |A 2 /A 3 | for solution points with tan/3 
3-15. 
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FIG. 7: Same figure as above but with tan (3 = 30 — 40. 
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FIG. 8: The correlation between solar neutrino mixing angle and |Ai/A 2 | for solution points with (a) tan/3 = 3 — 15 and (b) 
tan/3 = 30 - 40. 
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Stau LSP 






c (Ttn I 1 1 \ — 1 

b h n \P) — 1 






A = 39.5GeV 


m = 105.0 GeV 


Mi/2 = 671.0 GeV 


K 


8.40 x 10~ 6 


-7.66 x 10" 5 


-6.79 x 10 -5 


A, 


8.74 x 1(T 5 


-7.44 x 10~ 5 





e< 


-9.90 x 10~ 7 


-3.05 x 10~ 6 


-3.73 x 10~ 6 




1.01 x 10~ 6 


-1.31 x 10" 5 


-1.19 x 10" 5 


BR 


e 


/•* 


T 




28.9 % 


21.0 % 


50.0 % 


t 6 




~ 0.07 % 






m fl =278.6 GeV 


r = 


2.92 x 10~ 7 GeV 






~ 1.45 xlO~ 2 (Pb), 


= 1 TeV 



(Am§i, AmlO = (2.5 x 10~ 3 , 1.1 x 10~ 4 ) eV 2 
(sin 2 20 atm , sm 2 26 3ol , sin 2 26 chooz ) = (0.98, 0.77, 0.03) 
The decay length L ~ 6.74 x 10~ 8 cm 



TABLE I: A solution set with the stau LSP. The values of the trilinear couplings A; and Aj are taken as input parameters 
defined at the weak scale. 





Stau LSP 




tan/3 = 38.0 


sgn(/u) = -1 


H= -460.1 GeV 




A = 679.0GeV 


m = 233.4 GeV 


Mi /2 = 462.1 GeV 


a; 


1.53 x 10~ 6 


-1.48 x 10" 5 


-5.84 x 10" 6 


A, 


3.65 x 10~ 6 


5.41 x 10~ 6 





ti 


-7.76 x 10~ 6 


-2.81 x 10~ 5 


-2.30 x 10~ 6 


m 


8.17 x 10~ 6 


-3.38 x 10" 5 


-1.73 x 10" 5 


BR 


e 


A* 


T 


Uv 


11.5 % 


25.5 % 


47.2 % 


ib 




~ 16.0 % 






m fl =188.7 GeV 


r = 


7.99 x 10" 10 GeV 






~ 1.96 xlO~ 2 (Pb), 


Vs = 1 TeV 



(Am|i, Amli) = (2.5 x 10~ 3 , 4.6 x 10 _s ) eV 2 
(sin 2 20 atm , sin 2 2(9 soi , sin 2 26 chooz ) = (0.99, 0.74, 0.005) 
The decay length L ~ 2.47 x 10" s cm 



TABLE II: Same as the previous table but with large tan/3. 





Neutralino LSP 




tan/3 = 4.9 

-r±o — oo.yvjc V 


sgn(^) = -1 
rnn — 7 HpV 

//fcO — OOO. 1 Vjrc V 


At = -200.5 GeV 


A' 
K 

m 


-9.33 x 10" 9 
-5.63 x 1(T 5 
-1.21 x 1(T 6 
-2.01 x 10" 6 


-7.81 x 10" 5 
-7.35 x 1(T 5 
2.75 x 1(T 7 
-1.17 x 10" 5 


-7.56 x 10" 5 


1.83 x 10~ 6 
-8.74 x 10" 6 


BR 


e 


A* 


T 


vjj 
If 33 
ulfr* 


3.88 xl(T 2 % 
9.8 % 


47.0 % 
2.00 xl0~ 3 % 
16.6 % 


8.87 xl0~ 2 % 
26.4 % 




m.o=59.4 GeV 


r = 


7.14 x 10" 15 GeV 






.o.o ~ 4.90 xl0~ 2 (Pb), 


v/i = 1 TeV 



(Ami!, Amli) = (2.5 x 10~ 3 , 9.6 x 10" 5 ) eV 2 
(sin 2 20 atm , sin 2 2(9 soi , sin 2 26 chooz ) = (0.98, 0.99, 0.008) 
The decay length L ~ 2.76cm 



TABLE III: A solution set with the neutralino LSP which is lighter than the W boson. 





Neutralino LSP 




tan/3 = 4.2 


sgn(/i) = -1 


At = -442.3 GeV 




A = 98.5GeV 


m = 235.9 GeV M 1/2 = 380.9 GeV 


K 


-1.18 x 10 -8 


-1.44 x 10~ 4 


-1.47 x 10~ 4 




-8.57 x 10~ 5 


-9.71 x 10" 6 







3.24 x 10" 7 


-8.37 x 10~ 7 


-1.17 x 10" 6 


»7i 


-6.72 x 10" 7 


-1.73 x 10" 5 


-1.68 x 10" 5 


BR 


e 


A* 


T 


"33 




30.1% 




If 33 


3.14 x 10~ 2 % 


2.10 x 10 _1 % 


4.08 x 10 _1 % 




15.1 % 


19.4 % 


34.4 % 




2.35 x 10~ 2 % 


1.57 x 10 _1 % 


3.04 x 10 _1 % 




m.o= 165.8 GeV 


r = 


5.36 x 10~ 12 GeV 






~ 0.12 (Pb), 


^ = 1 TeV 



(Amii, Amli) = (2.5 x 10~ 3 , 6.5 x 10~ 5 ) eV 2 

(sin 2 26Un, sin 2 2(9 soi , sin 2 26 chooz ) = (0.99, 0.96, 0.001) 

The decay length L ~ 3.67 x 10~ 3 cm 



TABLE IV: A solution set with the neutralino LSP which is heavier than the W boson but lighter than the top quark 
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Neutralino LSP 






sgn(ju) = -1 


At = -832.3 GeV 




A —AS OripV 


m = 255.4 GeV M 1/2 = 715.3 GeV 


A' 


uo v i n — 6 

O.^O A _LU 


-4.99 x 1(T 5 


-7.26 x 10~ 5 


A, 


— ^ on v i n -5 


-6.22 x 1(T 5 





S 1 


2.12 x 1(T 6 


7.66 x 1(T 6 


8.53 x 1CT 6 


Vi 


-3.41 x 10" 7 


-5.84 x 10" 6 


-9.05 x 10 -6 


BR 


e 


/•* 


r 






11.6% 




If 33 


1.82X10- 1 % 


2.37 % 


2.95 % 


itW 


5.81 x 1(T 2 % 


1.65% 


3.13% 




14.8 % 


22.9 % 


37.0 % 




1.35 x 10 _1 % 


1.77 % 


2.19 % 




m -o=315.5 GeV 


r = 


2.14 x 10~ n GeV 






~ 0.14 (Pb), 


= 1 TeV 



(Am§i, Amli) = (2.5 x 10~ 3 , 2.2 x 10" 5 ) eV 2 
(sin 2 2(9 atm , sin 2 2(9 so! , sin 2 26 chooz ) = (0.91, 0.97, 0.05) 
The decay length L ~ 9.23 x 10 _4 cm 



TABLE V: A solution set with the neutralino LSP which is heavier than the top quark but lighter than two top quarks. 





Neutralino LSP 




tan/3 = 9.7 


sgn(At) = -1 


At = -1045.8 GeV 




A = 526.5GeV 


m = 308.5 GeV 


Mi/2 = 947.6 GeV 


a: 


9.00 x 10~ fi 


6.54 x 10~ 5 


-7.61 x 10" 5 


A, 


6.40 x 10~ 5 


7.51 x 10~ 5 





ti 


-3.49 x 10~ 6 


6.97 x 10~ 6 


-3.82 x 10~ 6 


Vi 


-3.56 x 10" 6 


1.96 x 10~ 5 


-2.12 x 10" 5 


BR 


e 


A* 


T 


vti 




97.8% 




1*33 


4.6xl0~ 3 % 


1.8xl0~ 2 % 


5.5xl0~ 3 % 


Ifitb) 


1.6 x 10" 3 % 


3.3 x 10~ 2 % 


4.3 x 10~ 2 % 




4.3 x 10 _1 % 


5.9 x 10 _1 % 


1.02 % 




3.39 x 10~ 3 % 


1.35 x 10~ 2 % 


4.05 x 10~ 3 % 




rn^o =417.0 GeV 


r = 


2.04 x 10" 9 GeV 






~ 1.04 x 10 _1 (Pb), 


■s/s = 1 TeV 



(Am|x, Arn^) = (3.6 x 10~ 3 , 3.2 x 10~ 5 ) eV 2 

(sin 2 26» atm , sin 2 26» so; , sin 2 28 chooz ) = (0.96, 0.98, 0.03) 

The decay length L ~ 9.66 x 10~ 6 cm 



TABLE VI: A solution set with the neutralino LSP which is heavier than two top quarks. 



